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ABSTRACT 

Context. With the discovery of the accelerated expansion of the universe, different observational probes have been proposed to 
investigate the presence of dark energy, including possible modifications to the gravitation laws by accurately measuring the expansion 
of the Universe and the growth of structures. We need to optimize the return from future dark energy surveys to obtain the best results 
from these probes. 

Aims. A high precision weak-lensing analysis requires not an only accurate measurement of galaxy shapes but also a precise and 
unbiased measurement of galaxy redshifts. The survey strategy has to be defined following both the photometric redshift and shape 
measurement accuracy. 

Methods. We define the key properties of the weak-lensing instrument and compute the effective PSF and the overall throughput and 
sensitivities. We then investigate the impact of the pixel scale on the sampling of the effective PSF, and place upper limits on the pixel 
scale. We then define the survey strategy computing the survey area including in particular both the Galactic absorption and Zodiacal 
light variation accross the sky. Using the Le Phare photometric redshift code and realistic galaxy mock catalog, we investigate the 
properties of different filter-sets and the importance of the u-band photometry quality to optimize the photometric redshift and the 
dark energy figure of merit (FoM). 

Results. Using the predicted photometric redshift quality, simple shape measurement requirements, and a proper sky model, we 
explore what could be an optimal weak-lensing dark energy mission based on FoM calculation. We find that we can derive the most 
accurate the photometric redshifts for the bulk of the faint galaxy population when filters have a resolution K ~ 3.2. We show that 
an optimal mission would survey the sky through eight filters using two cameras (visible and near infrared). Assuming a five-year 
mission duration, a mirror size of 1.5m and a 0.5deg^ FOV with a visible pixel scale of 0.15", we found that a homogeneous survey 
reaching a survey population of I^b=25.6 (lOcr) with a sky coverage of ~ 1 lOOOdeg- maximizes the weak lensing FoM. The effective 
number density of galaxies used for WL is then ~45gal/arcmin^, which is at least a factor of two higher than ground-based surveys. 
Conclusions. This study demonstrates that a full account of the observational strategy is required to properly optimize the instrument 
parameters and maximize the FoM of the future weak-lensing space dark energy mission. 

Key words. Photometric Redshift - Weak Lensing Surveys ~ Dark Energy - Cosmology 



1. Introduction 

With the measurement of the accelerated ex pansion of 
the U n iverse using Ty pe la S upernovae (Riess et al. 
(119981): IWood-Vasev et all (l2007l) : iKowalski et al. (2008) : 
iPerlmutter et al. (Il999h ). together with the flatness of the 
metrics derived from many CMB balloon-borne and sp ace ex- 

Sieriments (WMAP-7 years: Spergel et al. (2003); Komats u"et alJ 
20091) ). cosmology has entered a new era of precision mea- 
surements. The concordance Lambda cold dark matter model 
of the CMB and SNIa probes is also consistent with other 
probes (Baryonic Accoustic Oscillation (hereafter BAO) 



eg. lEisenstein eTaP (l2005h: iPercival et al.1 (l2010l) . cluster 
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counts eg. Takada & Bridle (2007), weak-lensing (hereafter 
WL) eg. Fu etjiL (2008)). This successful model has, however, 
reintroduced Einstein's controversial cosmological constant, 
which remains a mystery for fundamental physics. The contri- 
bution of the cosmological constant could be similar to that of a 
"dark energy" (hereafter DE) that would explain the observation 
of an accelerating Universe. Other theoretical models propose 
a change in the laws of gravity instead of adding an unknown 
"Dark Energy" component. Discriminating between the several 
DE solutions (kinder 2008) is the challenge of observational 
cosmology over the next decade. It has in particular motivated 
the preparation of future space-based missions such as JDEM, 
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the Joint Dark Energy MissiorQ on the US side ( for which 3 
concepts were in competition: SNAP DESTINY: i Morse et"ar 
(2004)) and on the European side the EUCLID mission, 
which represents the "merging " of the DUN^ and the SPAC: 
concepts. 

To go beyond our current limited observations of the 
Universe, we critically need new experiments that will pro- 
vide new and numerous observations of galaxies in the 
Universe to address the fundamental questions of cosmol- 
ogy. Different cosmology probes have been proposed to mea- 
sure the DE equation of s tate. These include in p articu- 
lar SN I a dbawson et al.ll2009l) . WL tomograph y (Mass ev et al l 



Hul 119991). anT3D-W L dKitching et al . 2007; Heavend 
Heavens et all l2006h. BAO dPadma nabhan & White 



200 9*). cluster counts (Marian et al 12009^), cluster strong 
lensing (Ju lio & Kneib 2009), and Alcock-Pazsinsky test 
dMarinoni & B uzzi 2010). The best approach is most likely to 
combine different probes, allowing us to minimize possible sys- 
tematic effects. 

WL has emerged as one o f the most effective cosmological 
probes ( Albrecht et al.l (|2006|) . see also the more recent JDEM 
FoM working group results dXlbrecht et al.ll2009l) ) as it is sensi- 
tive to both the geometry (through its dependence on angular- 
diameter distance ratio) and the growth of structure. The ob- 
served shape of a distant galaxy depends on the amount of mass 
distributed along the line of sight. To obtain the highest quality 
cosmological constraints, it is critical to derive accurate redshift 
measurements of all the galaxies for which one can measure their 
shape ( Massev et al. (20()7a)). In other words, any future WL 
imaging survey must address the question of the complementary 
redshift survey. We are presently unable to measure the galaxy 
redshifts of all the galaxies used in the shear estimation using 
spectroscopic technique. The only solution is to use photometric 
redshift. Although photometric redshifts have now been used for 
many years, the technique has mainly been developed using data 
available at various telescopes. However, very rarely has an in- 
strument or a survey been designed to optimize the photometric 
redshift measurement needed to reach a specific goal. 

Previous work aimed particularly at opti mizing photo- 
metric redsh ifts for future survey s include e.g., iBemtez et al.l 
(l2009i) and iDahlen et af] (|2008|) . which consider the filter 
properties, t heir nu mber and the photometry efficiency, and 
also Bordolo i et all (120 09). Schup (l2009h . lOuadri & Williamsl 
(12009). and Sheth & Rossii (|2010 C which evaluate the possible 
improvement of the photometric redshift technique using re- 
spectively, some work on likelihood functions, cross-correlation 
methods, close galaxy pairs, convolution, and deconvolution 
methods from a subsample of spectroscopic redshifts. A detailed 
study of the impact of photometric redshift errors on dark energy 
constraints was performed by Hearin et al. (2010) who gene ral- 
ized and extended the work of Bernstein & Huterej (lioiOh . It 
studies in detail the different types of photoz errors, their impact 
on dark energy parameters and the tolerances that will be use- 
ful in future survey design. The present paper extends the earlier 
work of Dahlen et al.. (2008,) and places the photometric redshift 
determination in the global context of the DE mission optimiza- 
tion. 



http://jdem.gsfc.nasa.gov/ 
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As we prepare future cosmological surveys, it is important 
to develop the optimal observational strategy and the photomet- 
ric data of a WL survey to maxi mize the prim e science of the 
DE mission based on the DETF (lAlbrecht et al. 2006) figure of 
merit. To achieve this goal, we use mock catalogs with realistic 
galaxy distributions as described in iJouvel et al.l (l2009l) (here- 
after Paper I) that is specifically designed to address this prob- 
lem. 

This paper is organized as follows. In section |2] we quickly 
summarize current photometric redshift techniques and charac- 
terize the likely photometric uncertainties of future WL mis- 
sions. We develop the WL requirements for future space DE mis- 
sions in section|3] In section|4] we investigate different filter con- 
figurations and underline the key characteristics of favored con- 
figurations. Section |5] investigates the impact of the blue-band 
photometry efficiency to help decrease the catastrophic redshift 
rate. 

Finally, in section |6] we explore the survey strategy in terms 
of a DE figure-of-merit (FoM) by investigating how the survey 
efficiency depends on the number of filters, the area of the sky 
surveyed, and the total exposure time per pointing. We discuss 
the results and possible improvements in section |7] 

Throughout this paper, we assume a flat Lambda-CDM cos- 
mology and use the AB magnitude system. 

2. Photometric Redshifts and Photometric Noise 

The photomeric redshift technique is to some extent similar to 
very low resolution (typically ~ 5) spectroscopy, but in- 
stead of identifying emission or absorption lines, it relies on 
the continuum of spectra and the detection of broad spectral 
features generally strong enough to be detected in visible and 
NIR filters. These features include "breaks" or "bumps" in the 
galaxy spectral energy distrib ution (SED) (Sawicki et al. 199^ 
iBolzonella et al.l 120001: iBenftezI 12000). Depending on the filter 
resolution, any spectral features that produce a change in colors 
can help the photometric redshift (hereafter photoz) determina- 
tion. 

There are three or four main spectral features that are par- 
ticularly helpful to the photoz procedure of which the most fun- 
damental are the Balmer break at ~ 3700A and the D4000A 
break. Additional useful characteristics are the Lyman break at 
9 1 2 A and the Lyman forest created by absorbers along the line of 
sight. However, the Lyman break only enters to the U-band filter 
at z ^ 2.5 and therefore only helps in breaking the color-redshift 
degeneracies for high redshift galaxies. In contrast the 1.6/im 
bump (Sawicki 2002) might be capable of breaking the color- 
redshift degeneracies of low redshift galaxies if a filter with cov- 
erage redder than 1.6jum is added to the filter set. 



2.1. Photometric redshift techniques 

There are two main types of methods that have been used 
to derive redshifts based on the photometry of objects: (1.) 
Empirical methods such as neural network (NN) techniques 
(Colhster & Laha vii2004j: I VanzeUa et al.ll200 4) and (2.) template 
fit ting methods su ch as the BPZ Bayesian phot o metric redshift 
of lBemtezI (|2000|). HvperZ of Bolzonella et al.l (l2000h . and Le 
PhareFlused in lllbert et aD(l2006i.i2009 ). Both methods includes 
two steps. The first step is the most critical in ensuring the ro- 
bustness of the photometric redshift estimate. For the NN tech- 
nique, this step is crucial. It uses a training set of galaxies from 
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which the NN learns the relation between photometry and red- 
shift. For the template fitting method, this corresponds to the 
calibration of the library of galaxy templates thereafter used in 
the redshift estimation. The template fitting method can work 
without this first step but it may then introduce some bias if the 
templates used are not representative of the galaxies for which 
the photometric redshift are measured. However, we aim to ob- 
tain unbiased photometric redshift measurements for many faint 
galaxies, it is essential that we calibrate the library of galaxy 
templates. Indeed, the calibration sample or the training set 
needs to be representative of the galaxy population for which 
we wish to find a redshift. The second step in both methods is 
the photometric redshift computation of the full galaxy sample 
from the photometry. The NN uses the complex function learned 
from the training set, while the template fitting method uses the 
calibrated library with a minimisation procedure to derive a red- 
shift estimation for each galaxy in a photometric catalogue. 

In this investigation, we use the Le Phare photometric red- 
shift code, which is based on the template fitting method. The 
code is applied to galaxies in a mock galaxy catalog that we 
describe in the next subsection. For each galaxy, the code de- 
rives a photometric redshift and a best-fit galaxy template using 
ax^ minimisation defined as 

n 

xlodel-Y^iKts-aFLdelV^'f (1) 
1=1 

where F'^^^^ and F'^^^^^^ are the observed and the template model 
fluxes inside a filter ; and cr' is the photometric error for this 
filter in a given survey configuration (as defined in section l231 l. 
Photometric errors play the role of a weight in the minimi- 
sation method and a is a normalisation factor The photometric 
redshift and best-fit template correspond to the minimum value 
of the;^:-^ distribution for a given simulated galaxy. 



2.2. CMC mock catalogue 

In Paper I, we developed realistic spectro-photometric mock 
galaxy catalogs. In this paper, we use one of those catalogs, 
flie COSMOS mock catalog (hereafter C MC), which was built 
from the observe d COSMOS data set jScoville et all 120071; 
ICapak etal.1 12008). This catalog use s the photometric redshift 
and best-fit template distribution of lllbert et al.l (l2009h . Using 
these two pieces of information, we calculate the theoretical 
fluxes of each galaxy in each band of a given survey configu- 
ration. We then draw an observed flux from a Gaussian distribu- 
tion based on the error estimate to simulate the observed galaxy 
photometric properties. The errors depend on the survey con- 
figuration, and the method used to calculate them is described 
in Section 12.31 We note that the mock galaxy catalog is pro- 
duced using the same set of templates utilised by the photometric 
redshift code. However, the representativeness of the calibration 
sample in the template fitting procedure is not the aim of thi s 
paper but will be studied in a future paper I Jouvel et al.l (1201 Oah . 
Thus we assume a "perfect" calibration in using the same library 
of templates for the development of the mock catalog and in the 
X^ procedure. Despite this being a very optimistic case, it pro- 
vides predictions and some results in the "optimal" case. 

The CMC assigns several emission lines to all galaxies in 
the catalog based on their [OH] fluxes, using the calibration of 
[Kennicutt ( 1998). The emission line fluxes are added to the flux 
derived from the continuum of each galaxy in the mock cata- 
logue. This creates a natural dispersion in the simulated magni- 



tudes, reflecting what will be observed in future real data. The 
bias that the emission lines will produce in the photometric red- 
shift estimate is one of the justifications for a photometric red- 
shift calibration survey (PZCS) ideally covering the same range 
in magnitude and redshift as the photometric galaxy catalogue. 
A wide and deep PZCS will help us to decrease the bias and dis- 
persion of the photometric redshift distribution using template 
calibration techniques. In optimizing the library of templates 
used in the photometric redshift analysis, we will be able to re- 
produce more accurately the diversity of the observed galaxy 
population including the impact of the emission line fluxes as 
shown in lllbert et al.l (12009,) . who found that their results are 
greatly improved where a spectroscopic galaxy sample is avail- 
able. The new version of the Le Phare code includes t he emission 
line fluxes in the library of templates as described in lllbert et al.l 
(2009). This last feature was a maj or impact in helpin g to im- 
prove photometric redshift results of lllbert et alj (l2009h . 

2.3. Typical noise properties for a space based survey 

In Paper I, we did not discuss in detail the typical photometric 
uncertainty caused by the instrument design and survey strategy. 
Since we wish to investigate the photometric redshift quality of 
future surveys, we now need to produce a realistic noise distri- 
bution for each galaxy in our catalog. To achieve this, we assign 
a photometric noise to each band that depends on the galaxy size 
and flux. Since we use electronic devices, the photometric signal 
is physically stored as electrons. Thus we express our formu- 
lae in terms of the number of electrons, which is proportional 
to the number of photons. We define esig„ai as the number of 
electrons produced by the galaxy flux. The photon noise can be 
described by a Poissonian statistic. Other sources of uncertainty 
originate in the instrument electronic devices and other astro- 
physical sources photons detected at the telescope. These stud- 
ies are space oriented so the main source of background noise 
comes from the Zodiacal light e^kv which is true in particular for 
a mission orbiting L2. The thermal radiation of the detector re- 
sults in a "dark current" ejark, while the reading of the detectors 
results in a read-out noise eRON described with a Gaussian statis- 
tic. We go through each of these four terms contributing to the 
noise in the Appendix. 

The signal-to-noise ratio including all the noise contributions 
is defined by 

S/N = 

^signal 

(2) 

where Nexpo is the number of exposures, T„bs the exposure 
time, and Npix the number of pixels taken in the flux error cal- 
culation. We took the RON to be e]^^^ - 6e^/pix and the 
dark current e^"^^ = Q.Q3e^ /pix/s for the visible detectors and 
^RON ~ 5e^/pix and e^"^.^ = Q.QSe^ / pix/ s for the NIR detec- 
tors. All parameter values are listed in the Appendix of table 
lA.ll These performances are achieved or expected in the near 
future from detectors of future DE missions. Thus, for each 
galaxy in each band, we calculate a S/N from equation |2] and 
compute an observed flux f^j^j' that includes a random noise 
drawn from a Gaussian distribution whose characteristics are 
(ju, cr) - {f^''li°,S /N), where is the noiseless or theoreti- 
cal flux value given by the CMC mock catalog. Thereby, using 
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the mock catalogs of I Jouvel et alJ (l2009t) and characteristics of 
future surveys, we compute realistic mock galaxy catalogs for 
future WL DE surveys including a redshift Zs, template model, 
galaxy fluxes, and uncertainties in each photometric band, in ad- 
dition to a galaxy size. More details about the calculation of 
the S/N are given in the Appendix. Following this noise pre- 
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Fig. 1. Iab magnitude as a function of exposure time for mir- 
ror diameters of 1.2m, 1.5m, and 1.8m using a pixel size of re- 
spectively 0.19",0.15", and 0.12" and a filter resolution of 3.2. 
The magnitude is calculated for four exposures (Nexpo = 4) 
assuming an exposure time by exposure written on the axis, 
a RON of - 6e^/pix and a dark current of e^"^^ = 

Q.Q5e^ / pix/ s. All parameters values are listed in the Appendix 
table lA.ll Magnitudes are calculated inside a circular aperture of 
1.4XFWHM. 



scription. Figure [T] shows the I-band magnitude as a function of 
exposure time for a given S/N ~ 5 and 10 in blue and green, 
respectively, and for mirror sizes of 1 .2m (small-dashed lines), 
1.5m (large-dashed lines), and 1.8m (solid lines). These values 
are derived assuming an obstructed telescope design with a mir- 
ror size for the secondary of 60% of the primary mirror. This 
shows for example that a 1.5m telescope and a survey strategy 
of four exposures of 200s (800s of total integration time) reaches 
a magnitude of /ab=25.8 {S /N ~ 10) for a galaxy source of 
FWHM"'", = 0.20, where FWHM"'"; is the observed FWHM 

gal gal 

of a galaxy. Magnitudes are computed inside a circular aper- 
ture of 1.4xFWHM. The stars in gold represents the exposure 
time needed to reach the COSMOS completeness for different 
telescope diameters calculated using our noise prescription. The 
magnitude calculation is described in the Appendix. 

To obtain an accurate WL measurement, it is safe to use the 
galaxies whose FWHM are larger than 1.25x[FWHM(ePSF)] 
and > 10, where the ePSF is the effective PSF of the tele- 
scope defined in section 13.21 The COSMOS WL analysis used 
a criterion of 1.6x[FWHM(ePSF)] and a S/N>10, but we hope 
that an image analysis technique of higher quality will improve 
the COSMOS limit in the future. The choice of a factor of 1.25 



although an arbitrary criterion, allows us to easily compare dif- 
ferent survey designs by using a simple size cut as a quality cut. 

Figure |2] shows the number density of galaxies reached us- 
ing these criteria for a primary mirror size of 1.2m, 1.5m, and 
1.8m. In decreasing the primary mirror diameter by 0.3m, the 
galaxy number density is reduced by 13gal/arcmin^ when go- 
ing from 1.8m to 1.5m, and 18gal/arcmin^ when going from 
1.5m to 1.2m. We choose to use a pixel scale that varies with 
the mirror size to ensure an equal sampling of the effective PSF. 
We choose, respectively, a pixel scale of 0.19" for a mirror size 
of 1.2m, and 0.15" for 1.5m and 0.12" for 1.8m. Figure |2] also 
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Fig. 2. Effective number of galaxies as a function of exposure 
time for a mirror diameter of 1 .2m, 1 .5m , and 1 .8m using a pixel 
size of respectively 0. 1 9",0. 1 5", and 0.12" and a filter resolution 
of 3.2. All parameters used to produce this figure are listed in the 
Appendix Table lA.fl 

shows that the quality cut based on galaxy size produces a loss 
of 3, 5, and 6gal/arcmin^, respectively for a primary mirror di- 
ameter of 1.8m ,1.5m, and 1.2m. We note that the loss is more 
significant for smaller mirror sizes. This is due to the relation- 
ship between the mirror size and the pixel scale. Smaller miiTor 
have larger pixel scales, which makes the quality cut on galaxy 
size more stringent. However, this is a small loss compared to 
the 31gal/arcmin^ that one loses when going from a mirror size 
of 1.8m to 1.2m based on 4 exposures of 200s. Using the expo- 
sure time needed to reach the COSMOS completeness (shown in 
Figure|2]i, we have a galaxy density of 71gal/arcmin2. This de- 
fines an exposure time-density domain in which the COSMOS 
catalog and the CMC are complete. The dashed gold region cor- 
responds to areas where the CMC catalogues produced are in- 
complete. In these areas, conclusions may be affected by the in- 
completeness of the CMC. 

2.4. How fo characterize the photometric-redshift quality 

We call the redshift coming from the input mock catalog the 
spectroscopic redshift Zs, while the photometric redshift Zp cor- 
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responds to the redshift calculated by the photometric redshift 
technique. Considering a photometric redshift distribution, we 
define the "core of the distribution" as the galaxies for which 
\zp - zs\ < 0.3 and the "catastrophic redshift" as the galaxies 
outside the core. We did not include the division by 1 + zi since 
we had not intended to produce results to be used in WL anal- 
yses, but to instead assess the photometric redshift quality. We 
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Fig. 3. Illustrative diagram of photometric versus spectroscopic 
redshift, where we identify the quantities assessing the photo- 
metric redshift quality. 



define some characteristic numbers that we use to quantify the 
quality of a photometric redshift distribution: 

- cTcore, the dispersion of the core distribution defined as 

cr(\Zp-Zs\<03). 



2AA . A confusion between the Lyman and Balmer breaks 

This break confusion is represented by the four purple lines. In 
an ideal case, if the photometric redshift procedure identifies a 
higly accurate redshift zp - zs, we find that 

^break-obs - (1 + Zs)Ai,reak-rf , (3) 

where Ahreak-obs is the observed wavelength of one of the breaks 
used in the photometric redshift procedure, Abreak-rf is the rest- 
frame wavelength of the same break, and zs the spectroscopic 
redshift of the galaxy. In the case of a catastrophic redshift zp - 
Zcata, this equation becomes 



Abreak-obs — (1 Zcata) Abreak-cai 

We can then write: 

/I , Abreak-rf ^ 
■Zcata = (1 +ZS)- 1 



(4) 



(5) 



(break — rf, break — cata) — 



(6) 



Abreak-cala 

We define four line couples (break - rf, break - cata) that are 
sources of the color degeneracy producing the catastrophic red- 
shifts, where break - rf is the real feature and break - cata is 
the wrong feature found by a photoz code : 

Ly a Balmer \ 

Ly a D4000 

Ly break Balmer 

Ly break D4000 , 

These couples used in equation |5] define the four purple lines in 
the zp-zs plane where the catastrophic redshift happens with the 
highest probability. 

This confusion occurs for both low and high redshift galax- 
ies, generally at z < 0.5 and z > 2.5, depending on the wave- 
length range available to the instrument. A wide wavelength 
range going from U to K band would avoid most catastrophic 
redshifts by using both the U-band and NIR photometry. The 
Balmer break can be followed at all redshifts from the V-band 
photometry (z ~ Balmer break~ 4000 A) to H-band photome- 
try (z ~ 3 Balmer break~ 16000 A). However, it can be misiden- 
tified as the Lyman break leading to the creation of catastrophic 
redshifts. This misidentification can be avoided by using deep 
U-band photometry. 

The break confusion will generally produce a double peak 
in the redshift probability distribution of low-redshift galaxies 
< Zi < 0.5, one at the correct redshift, and one at higher red- 



distribution defined as i-i(\Zp - Zs\ < 0.3). 
ricore, the number density of galaxies inside the core distribu- 
tion. 

n"'"", the number density of galaxies with a photoz of high 
confidence that we defined as ^^^'^'z < 0.5. 
ncata, the number density of galaxies with a catastrophic red- 
shift defined as \zp - Zs\ > 0.3. 

n'™^', the number density of galaxies with a photoz of high 
confidence being catastrophic redshifts. 



- //(-ore, the bias measured from the mean or median of the core shift 3.5 < Zp < 4, which corresponds to a "catastrophic red- 
shift". Hence, the derived photometric redshift distribution can 
be biased having an excess of galaxies with 3.5 < z„ < 4, which 
will strongly perturb the DE parameter estimation (iHuterer et al.l 
2006). In section |5] we investigate more quantitatively the gain 
of an eflicient U-band in minimizing the break confusion. 

2.4.2. An inaccurate template fitting 

The photometric redshift dispersion and biases depend on the 
quality of the photometry of galaxies (which can be affected 
by instrumental defects or crowded fields). Deeper photometry 
helps to provide higher accuracy photometric redshifts at a given 
magnitude. The galaxy color accuracy is higher with deeper pho- 
tometry and the weight of the fit given by the S/N is higher, 
which both decrease the dispersion and possible biases in the 
photometric redshift estimate. In addition, a slight filter calibra- 
tion error is enough to bias the photometric redshift distribution. 
A way in minimizing the dispersion and biases of the photomet- 
ric redshift estimate is to optimize the resolution of the photo- 
metric bands. We explore this solution in section l4~2l 



Figure |3] is an illustrative density diagram showing pho- 
tometric versus (vs.) spectroscopic redshift with the core of 
the distribution being located inside the two red lines, and the 
catastrophic redshifts outside these red lines. Following this 
definition of catastrophic redshift, there are two kinds of "bad" 
redshift. The galaxies surrounding the red lines and the galaxies 
situated close to and within the two purple lines. There are two 
main reasons for the redshift procedure to fail, which we discuss 
below. 
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3. Weak lensing survey key parameters and 
definitions 

To reach the goal of precision cosmology, it is essential to opti- 
mize the instrument design and survey strategy, which both im- 
pact the quality of the WL results. The present section aims to in- 
troduce the quantities used in the DE parameter estimation such 
as: (1) the galaxy number density which is a function of the ex- 
posure time and the photometric redshift quality (2) the survey 
area including the impact of the Galactic absorption (3) the pixel 
size which impact the quality of the photometry and the shape 
measurement (4) the minimum exposure time to be in the photon 
noise regime. 

3.1. Weak lensing dark energy parameter list 

One of the possible way of constraining DE is the WL tomogra- 
phy de scribed in either Hu & Jain ( 2004) or Amara & Refregier 
(l2007h . This method divides the source distribution in redshift 
slices, thus requires that accurate and unbiased photometric red- 
shifts be available for most galaxies. A number of factors affect 
the FoM of this technique including, (1) the number of galaxies 
useful to the WL measurement, (2) the systematic errors in the 
shape measurement, and (3) the errors and biases in the photo- 
metric redshift distribution. 

The FoM formalism of the iCosmo package (iRefregier et alj 

I2OO8 ) is based on the WL tomography method. Using the galaxy 
densities defined from the photometric redshift results from our 
mock catalogs and the FoM from iCosmo, we look at the impact 
of the photometric redshift quality on the DE parameter esti- 
mation. We assume a flat cosmology where the fiducial values 
of cosmological parameters are (Q^, wo, Wa, /z, Q;,, erg, Qa) - 
[0.3, -0.95, 0, 0.7, 0.045, 0.8, 1 , 0.7]. We compute the FoM of the 
(wo, Wa) DE parameters in marginalising over the other cosmo- 
logical parameters and using five tomographic redshift bins that 
have been found to provide the most accurate FoM (Sun et al. 
l2()09l). The redshift di s tributi on f ollows a distribution de scribed 
in lSmail&DickinsonI (Il995h and lEfstathiou et all d 19901) 



N(z)^z"exp-(l.41z/z,nedf 



(7) 



with parameters a,fi = [2, 1.51 fo llowing the COSMOS red- 
shift distribution fit of Massev et al.l l2007b). The boundaries of 
the tomographic redshift bins are calculated to produce an equal 
repartition of the number of galaxies in each of the five redshift 
bins. To calculate the FoM, the key numbers that we derive from 
our mock catalogs are 

- Ngai, the galaxy number density of galaxies that satisfy 



<e S/N, 



I-band 



1- — — / - - ( — ii/itin — 1 
+ zp 

FWHM,,„, > 1.25 X FWHM.psF, 



(8) 



where ePSF is defined in section [X2l and e = (0.1,0.5) is a 
parameter that defines the quality of the photometric redshift. 

- Zmed is the median of the photometric redshift distribution of 

Ngal- 

- J?l is the survey area derived from the instrument field-of- 
view and the survey strategy, explained in section [341 

The number of objects Ngai depends on the photometric redshift 
error criteria that we assume, which are parametrized by e (stud- 
ied in section|6]l, the primary mirror size {D\), and the pixel scale 
(/?'"), which enter in the definition of the effective PSF (ePSF) 



discussed in section [X21 and in the photometric uncertainties de- 
scribed in section 12.31 In section 13.21 we define a maximal and 
an optimal pixel size by means of their impact on the size of the 
effective PSF, which determines the useful number of galaxies: 
Ngal. In section 13.31 we define the minimum exposure time in 
the photon noise regime depending on the instrument parame- 
ters. We then study in section [3741 the survey area J[ taking into 
account the Zodiacal light and Galactic absorption, which both 
depend on the sky position. 

3.2. ePSF: Effective PSF of the telescope and optimal pixel 
scale 

The future observation strategy is to survey a large fraction of the 
sky. This would be easier using large pixels typically of the order 
of the PSF size, which is a function of the mirror size (see Table 
lA.U : this would help to optimize the area versus observation 
time without under-sampling the PSF too much, which would 
affect the quality of the WL measurement. In this section, we 
define the pixel scale to be used in the calculation of the noise 
properties. 

3.2. L Formalism 



Using the formalism o f lHigh et al.l (I2OO7I) . the observed galaxy 
shape lobsiff) is expressed as the convolution of three components 
of the intensity profile of galaxies, the pixel response p{0) and 
the PSF of the telescope 



I"''' (9) = is'^'''-'y{Q) * PSF {6) * pie). 



(9) 



The theoretical PSF size PS F{0) corresponds to the size of the 
Airy disk. Its size is a function of the wavelength A and the mir- 
ror size on the basis of the relation 



PSF{e) = 2.44—. 



(10) 
(11) 



Similarly the full width half maximum of this PSF is defined by 



¥SNY{M{Airy disk] = l-02.-^, 



(12) 



where D\ is the diameter of the primary mirror. The PSF and 
pixel response introduce systematics in the WL measurement 
and need to be extracted from the galaxy shape before doing 
any lensing calculations. For this purpose, we use point-like 
sources such as stars to correct for both PSF circularization 
and anisotropic deformation. Thus, we define the effective PSF 
(ePSF) corresponding to the star intensity profile, which is the 
convolution of the PSF and the pixel response that will be ob- 
served on telescope images 



ePSF{0) = 5(6) * PSF{ff) * p{ff). 



(13) 



This ePSF corresponds to the resolution of the instrument or the 
smallest size resolved by the telescope. To obtain a rough esti- 
mate of the size of the ePS F, we assume Gaussian distributions 
for the pixel response, the PSF of the telescope, the jitter, and 
the pixel diffusion. The jitter and pixel diffusion also affect to 
the size o f the observed PSF and need to be taken into account 
(iMa et alji200&) . Thus we define the effective PSF expressed in 
arc sec as 



ePSF = yl{PSF(e)Y + 0.2 p^ 



(14) 
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where p - p"^/"' is the pixel scale (visible or IR camera), crj 
represents the jitter of the telescope and cr^ the diffusion of the 
pixel. The pixel diffusion varies as a function of the pixel size 
(p) and has a typical value of crj =0.04" for a pixel size of 0. 1", 
which is equal to a diffusion of 0.4 pixel (see Table lA.ll . 



3.2.2. Maximal and optimal pixel scale 

Extrapolating from the results of iHigh et al.l (l2007h . we define 
the maximal pixel scale. In the context of a DE WL survey, 
[High et al. (2007) defined an optimal pixel size for a primary 
mirror size of 2m to be 0.09" with one exposure at 0.8//m. This 
pixel scale slightly undersamples the ePSF. However, the pixel 
scale can be increased if a combination of sub-pixel dithered 
images are used. Different techniques can be used to com- 
bine the sub- sampled images such as the drizzling technique 
jFruchter & H ook 2002) working in real space or the method 
proposed by iLauer (.1999.) which works in Fourier space. To 
recover the loss of information caused by the ePSF under- 
sampling, the minimum number of exposure has to be N,„in - 

^priT j J where p"-^ is the under-sampled pixel scale and p'" the 

optimal pixel scale for one exposure. iHigh et al.l (l2007l) showed 
that for a primary mirror size of 2m, a pixel scale of 0.16" us- 
ing four perfectly interlaced images is a good alternative to one 
exposure with 0.09", if assuming a perfect image reconstruction 
from the four dithered exposures. In terms of PSF sampling, this 
would allow us to undersample the PSF by a factor ^ of 

PS F size (Di ^ 2m) 0.14 



MTF of ePSF 



MTF of PSF 



(15) 



Using the under-sampling factor ^ and a 1 .5m telescope, we find 
a pixel scale of p"^ ^ PSF size (Di - 1.5m)/f ^ 0.2", and for a 
1.2m telescope p"^ ^ 0.25". 

3.2.3. Pixel scale estimated with the MTF (modulation 
transfer function) 

Another way to define the optimal pixel scale of our configura- 
tion (Di = 1.5m and Nexpo = 4) is to apply the Nyquist-Shannon 
theorem. This theorem says that a function is completely deter- 
mined if it is sampled at ' where f,„ax is the highest fre- 
quency of the given function. We thus trace the Fourier trans- 
form coefficients of the ePSF and check that a pixel scale of p''" 
is given by p''" > 57" ■ Figure |4] shows the modulation transfer 
function MTF of the PSF and ePSF of a 2m mirror, where ePSF 
is defined in equation [14] The MTF corresponds to the coeffi- 
cients of the Fourier transform, that we trace as a function of the 
frequency. This shows that using the instrumental characteristics 
we defined, a pixel scale of 0.15" is a good choice, if we assume 
a perfect information recovery using the dithering technique. A 
pixel scale of 0.075" samples the ePSF well and allows us not 
to lose any information at any frequency. We translate this in- 
formation to a mirror of 1.5m using Figure[5| This figure shows 
the ePSF sampling in an I-band filter as a function of the primary 
mirror diameter for different pixel scales. The ePSF diameter is 
defined as two times the ePSF radius in equations [T2| and [T4| 
in section [3] This figure shows that a pixel scale of 0.15" for a 
primary mirror size of 2m is equivalent to a pixel scale of 0.2" 
for a primary mirror size of 1.5m in terms of ePSF sampling. 
We thus define 0.2" pixel as the maximal choice for future WL 
surveys, assuming a perfect information recovery with a perfect 
half-pixel dithering. With the perfect sub-pixel dithering, we can 
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Fig. 4. (right) MTF[PSF] and (left) MTF[ePSF] at ^ ^ 8000A, 
for a mirror size of 2m. The green line stands for a pixel scale of 
p"" = 0.0375", cyan p"' = 0.075", and blue p"' = 0.15". The 
ePSF comes from equation with a pixel scale of - 0.075", a 
jitter cTj ~ 0.04 and a diffusion crj - 0.04. The colors represents 
different pixel scales studied for the optical camera in arcsec. 
The pixel scale of the NIR camera is explained in equation [TSl 
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Fig. 5. ePSF sampling as a function of primary mirror diam- 
eter for different pixel scale at 800nm. The pink and orange 
points represent respectively the WFPC2, and ACS camera of 
the Hubble Space Telescope. The purple point represents a mir- 
ror size of 1.5m and a pixel scale of 0.2". 



hope to recover the information to reach O.F'/pixel with a sam- 
pling of 3 pixels for the full ePSF equivalent to 1 .5 pixel over the 
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FWHM[ePSF]. Such a c onfiguration will be discussed in more 
detail in a separate paper (I Jouvel et alJl2010bl) . 

The maximum pixel scale of 0.2" for a 1 .5m telescope should 
be considered as the upper limit while a safer solution, which we 
suggest is "optimal", has a pixel scale of 0.15" (two pixels sam- 
pling of ePSF), which is comparable in terms of the sampling of 
the ePSF to the sampling of the WFPC-2 camera of the Hubble 
Space Telescope. 

A higher samp l ing rate was proposed by 
iPaulin-Henriksson et al.l (l2008h to reach a good ePSF sam- 
pling and measure accurately this ePSF, allowing excellent WL 
measurement. However, enlarging the pixel size also provides a 
larger field of view (for a given number of detectors), which is a 
key parameter in the FoM determination. As we expect the ePSF 
of a s pace mission at L2 to be extremely stable (iBernstein et al.l 
I2009h thus well constrained, a full optimization of the pixel size 
that takes account of full observation strategy and the final WL 
FoM must be investigated before committing to a final design. 

3.2.4. Discussion 

Figure |6] shows the sampling of the full ePSF in the top figure 
and the FWHM of this ePSF as a function of wavelength using 
pixel scales of 0.05",0.07",0.15",0.2", and 0.25" respectively in 
blue, green, cyan, gold, magenta, and brown. These pixel scales 
correspond to the visible CCD detectors. For simplicity, we as- 
sume that the NIR detectors share the same focal plane so that 
the pixel ratio is just the ratio of the pixel physical size 



p"' - [ratio size] X - 



xp"" = 1.71 X 



(16) 



where p"' and are, respectively, the pixel scale of the NIR 
detector and the visible detector. We consider for the NIR detec- 
tors the physical size of 18;um and 10.5yL(m for visible detectors 
(LBNL CCD). 

We note that the ePSF is similarly sampled in the NIR wave- 
length range. The PSF size is proportional to the wavelength A 
such that PSF cc A/Di, where Di is the primary mirror size, that 
allows a higher sampling. However, a large PSF causes a sub- 
stantial decrease in the galaxy number density. The WL analysis 
makes use of the shape of galaxies. Thus, one has to make a cut 
in the galaxy size to use only the galaxies whose shapes are not 
contaminated by the instrumental PSF. As an illustrative exam- 
ple the PSF size at A - 1200nm with a mirror size of 1.5m is 
equivalent to that at A - 800nm for a mirror size of Im. Even 
if the count slope difFeres between J band and I band, it will de- 
crease the galaxy number density significantly, suggesting that 
the WL measurement should be conducted more efficiently in 
the visible bands. 
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Fig. 6. (Top) FWHM of the ePSF in arcsec and (Bottom) ePSF 
sampling as a function of wavelength in A for a mirror size of 
1.5m. 



3.3. Exposure time 

To establish the optimal DE FoM, we need to define a "mini- 
mal" exposure time for WL, which is a combination of three de- 
pendent factors: 1) the photon noise, which must dominate the 
detector noise for typical galaxy photometry and shape measure- 
ments; 2) the exposure time should be small enough to cover the 
largest possible area of available sky, but long enough to obtain 
a good photometric redshift distribution; 3) To reach a homo- 
geneous survey across the sky, the exposure time should be ad- 
justed depending on the Zodiacal light and Galactic absorption. 
This is particularly important at visible wavelengths where the 
WL measurement will be conducted. 



To study the impact of the filter set properties on the photo- 
metric redshift quality, we first have to define a minimal expo- 
sure time beyond which the photon noise dominates the detec- 
tor noise. Thus, for a given exposure, we extract the minimum 
exposure time T"™^ from which the read-out noise becomes sub- 
dominant in using the denominator of equation |2] 

^signal "t" ^ pixT fjjy^e f;i^yi sec — ^ pix^RON ^P^^'^ohs^^'^^^^ (1'^) 

where the left-hand term is the photon noise contribution from 
the Zodiacal light and a galaxy, respectively, e signal and ^5^, and 
the right-hand term is the detector noise with the read-out noise 
(RON) and the dark current (for more details, see the Appendix). 
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Fig. 7. Signal-to-noise ratio as a function of observation time for 
different I-band magnitudes. The telescope characteristics are 
listed in Table lA.fl We use the I-band filter for the blue, green, 
magenta, and red curves. The cyan curve uses the properties of 
the ACS camera in the F814W filter The square cyan dot on 
this last curve represents the COSMOS survey with an observa- 
tion time of 507s by exposure. The dotted gold line separates 
the RON (read-out noise) dominated regime (exposure time less 
than 135s) from the photon noise dominated regime at longer 
exposure times. The dotted cyan line represents the same as the 
dotted gold line for the ACS camera. 



Figure |7] shows the S/N as a function of the observing 
time at different I-band magnitudes of 24,25,26, and 27 (red, 
magenta, green, blue respectively) using a filter resolution of 
% - 3.2, a pixel size of 0.15", and a mirror diameter of 1.5m. 
The cyan curve corresponds to the ACS-COSMOS expectations 
for this noise simulation at a magnitude of 26 in the F814W 
filter with the ACS pixel scale of 0.05" and the HST mirror 
size of 2.4m. The square point represents the COSMOS survey 
(with four exposures of 507s) performed using the ACS camera, 
which has a RON of e^'^^ = le^ I pixel and a dark current of 
^^"r/t ~ \2>e^ I pixel I hour as stated in iKoekemoer et al] (l2007h . 
We find a limiting magnitude of AB(F814W)=26.6 (5.2o-) for a 
galaxy size of 0.21" effective radius. This estimate is very close 
to that of Leauthaud et alj (120071) . who find a limiting magnitude 
of AB(F814W)=26.6 (Scr) for a galaxy eff^ective radius of 0.2". 

We note a change of slope for the solid line curves denot- 
ing the detector and photon noise regime. For short exposures, 
the read-out noise dominates the denominator term and the S/N 
grows proportionally to the exposure time as shown by the ma- 
genta dot-little-dashed line. Thus, e signal °^ Tots and the S/N is 



SIN • 



^signal 



oh Si 



(18) 



2 

RON 



where e^oN - 6e^/pix is not time-dependent. For long expo- 
sures, the S/N varies as the square root of the observing time as 
shown by the magenta dot-long-dashed line 



S/N 



'^signal 



obs 



P Tab s 



(19) 



obs 



where a holds for the galaxy and sky photons, whose fluxes are 
a function of Tohs, and fi can be deduced from equation I A. 41 in 
the Appendix. 

Equation[l7]can be simplified to define a minimum exposure 
time at which the sky noise equals the collective detector noises 
(the dark current and read-out noise) 



" obi 



e 



sky I sec ~ 



'RON ^ ^ obs • 



dark-) 



(20) 



where e^ky - e sky/ sec Tohs- We find a minimum exposure time 
of 135 sec for a 1.5m mirror diameter with a 0.15" pixel scale. 
This number defines the minimum exposure time for which a 
WL survey is optimal. It is interesting to raise the S/N as long 
as we are in the detector noise regime and S /N cc Tabs- We also 
note that current software assumes that the noise properties fol- 
low Poisson statistics, which is true in the photon noise regime. 
However, in the detector noise regime, the noise properties fol- 
low Gaussian statistics. If not taken into account, this will impact 
the galaxy properties calculated from the software as well as the 
galaxy extraction. We return to the observation time in section|6] 
where we study its impact on DE parameter estimations. 
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Fig. 8. Minimum exposure time for WL surveys as a function of 
wavelength integrated in filters of resolution K = 3.2. The min- 
imal exposure time is defined in equation |20] and corresponds 
to the photon noise dominating detector noises. It takes into ac- 
count the mirror size, the pixel scale, and the filter efficiency. 
The thickness of the lines grows with the mirror size. 

Figure |8] shows the two dependences on both the mirror 
and pixel scales of the minimum exposure time as a function of 
wavelength integrated in filters of resolution 1^-3.2. This fig- 
ure shows that a smaller mirror and pixel scale requires a longer 
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minimum exposure time for the photon noise to dominate the 
detector noise. We note that a 1 .2m mirror diameter with 0. 1" by 
pixel requires a minimum of 700s exposure to be photon-noise- 
dominated. It is possible to decrease the minimum exposure time 
needed by using filters that are broader than our optimal resolu- 
tion of !R = 3.2. However, this would reduce the photoz quality 
(as shown in section |4]l and may jeopardize the PSF color cor- 
rection. The shape of curves reflect the logarithmic width of the 
filter set configuration (shown in Figure 12111 and the drop in the 
detector efliciency at blue wavelengths (as studied in FigurefTSll. 
This also explains the decrease in the sky background magnitude 
at blue wavelengths shown in Table [T] This noise magnitude is 
the Zodiacal light flux (explained in the Appendix) integrated 
within the photometric bands of the eight-filter set without tak- 
ing into account any instrument characteristic other than the fil- 
ter efficiency. In the table, 77 represents the whole transmission 
including filter transmission, mirror reflectivity, and detector ef- 
ficiency. 

We note that a fixed exposure time of 200s allows us to use 
optimally the information contained in almost all bands, except 
the two bluest bands (which would ideally require longer expo- 
sure times). For simplicity, we thus choose to use this exposure 
time to study the photometric redshift quality as a function of the 
resolution of filters in section |4] 



Table 1. Noise magnitude (Zodiacal light) in mag/arcsec^ for the 
eight-filter set configuration and the total telescope throughput 77 
in each band. 



To reach a homogeneous data quality, we need to adjust the 
exposure time of the survey as a function of the pointing posi- 
tion on the sky (a,6). We can define the exposure time factor 
needed to reach the intrinsic magnitude limit as a function of the 
coordinates (assuming that the survey is photon-noise-limited) 

zodi(a, 6) 



tf(a,6) - 



trans(a,6y-' 



(21) 



where tf is the exposure time required to reach the desired limit 
displayed in Figure [TOl zodi is the Zodiacal background light 
level plotted in Figure HI and trans is the Galaxy absorption 
defined by 

trans(a,6) = l0-"-4'^-'<'>'-«, (22) 

where is the extinction map at wavelength A due to Galactic 
dust from .Schlegel et al.. (,1998.) shown in Figure|9] 
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3.4. Galactic absorption, Zodiacai iigiit, and survey area 

For most cuiTent surveys (COSMOS, CFHT-LS, RCS2), the cor- 
rections for Galactic absorption and Zodiacal light are not diffi- 
cult to make. These surveys cover relatively small fields and are 
generally located at high Galactic latitudes. This will not be the 
case for the next generation of weak lensing surveys, which will 
be limited by both Galactic absorption and Zodiacal light varia- 
tion. 

In general, the overall number of galaxies grows faster when 
surveying wider fields rather than going deeper in smaller fields 
(which reflects the small gradient of the galaxy count slope). 
Future cosmological surveys should cover more than ten thou- 
sand square degrees as advocated by Amara & Refregier (2007), 
who demonstrated that DE constraints grow proportionally to 
the number of galaxies. However, when reaching such wide ar- 
eas, the impact of Galactic absorption and Zodiacal light vari- 
ation has to be accounted for in the survey strategy or it will 
otherwise severely affect the photometry quality, leading to 
degradation of the DE co nstraints. This was not addressed in 
lAmara & Refregied (l2007h . 



Fig. 9. (Top) Zod iacal background lev el as a function of eclip- 
tic latitude from iLeinert et al.l (l2002h . assuming the telescope 
viewing angle is between 70 and 110 degrees from the Sun. 
(Bottom) Galactic ab sorption as a function of wavelength from 
ISchlegelet al1([T998l) . 

In section |6l we use this model to investigate the fraction of 
the sky a telescope should survey to optimize the cosmological 
constraints. To do this, we need to define the telescope charac- 
teristics, such as the number of filters rif for each of the two 
cameras (ricam - 2: one infrared with HgCdTe detectors and one 
visible with CCDs - assumed here to have exactly the same field 
of view), the number of exposures Nexpo - 4 per filter, and the 
observation time Tots per exposure. These parameters define the 
survey configuration. 

Using the characteristics of the survey configuration, we de- 
fine the minimum exposure time for each pointing as 

"^obs/pting ~ Tohs/cam Ne.xpo n fleam, (23) 

where N^xpo is the number of exposures (see Table I A. fl . rtf/cam 
the number of filters by camera, and Tobs/caw the individual im- 
age exposure per filter for a given camera. 
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Exposure factor Survey duration in years 



Fig. 10. Cumulative distribution of the sky coverage for a given 
exposure time factor, i.e. sky coverage with exposure factor less 
than the given exposure factor, at several wavelengths. 



The exposure time for a given position on the sky (a, 6) is 
then defined as 



Tobs/ptingia, S) - tf(a, 6) x T[ 



obs/pting^ 



(24) 



where, for simplicity, the exposure time factor is computed at a 
wavelength of 800nm, which corresponds to the wavelength of 
the weak lensing measurement, and is applied globally. 

We define the survey area for a given camera field-of-view 
(FOV) and a total mission time Tmissmn as 



FOV X T„, 



1 '^obs/plingiP) 



(25) 



n) 



We note that Tmission includes a survey efficiency ( - 0.7, which 
accounts for observation overheads (telescope slewing, guide 
star acquisition, read out time, ...) and data transmission. The 
denominator is the mean observation time over the whole field 
surveyed and is calculated iteratively. Figure [11] shows a cumu- 
lative distribution of the sky coverage as a function of the total 
mission time of a survey in years. This figure uses a survey strat- 
egy as defined above consisting of changing the exposure time 
as a function of the Galactic absorption strength on the sky area 
observed. The black line includes the Galactic absorption and 
Zodacal light variation, while the blue line does not. 

This observation strategy ensures a uniform photometry 
quahty over the whole survey by spending more time on point- 
ings closer to the Galactic plane, as shown in Figure [121 This 
figure is a version of the sky map showing the exposure time re- 
quired to achieve a S/N of 10 for an Iab - 25.6 galaxy. The scale 
is in seconds. This assumes that four exposures of the time listed 
were taken. The time spent then depends on the field location on 
the sky and is determined with the exposure factor. 

We note that the survey area J[ scales as FOV x T mission, 
hence a reduction in the camera FOV thus reduction in the num- 
ber of detectors can be compensated for by a longer mission 
time. 



Fig. 11. Cumulative distribution of the sky coverage as a func- 
tion of the mission duration in years using the survey character- 
istics described in Table lA. fl at 8000 A. The black curve includes 
Galactic absorption and Zodiacal light variation, while the blue 
curve does not. 




■ 



■ 



Fig. 12. Exposure time map of the sky (in seconds) needed to 
reach S/N=10 for a Iab - 25.6 galaxy at a wavelength of 8000 
A(using the survey characteristics described in Table lA. ll . 



4. Filter resolution studies 

To optimize the survey strategy, we must carefully study all pa- 
rameters that affect the galaxy photometry. Using the noise prop- 
erties defined in section 1231 the optimal pixel scale, and the ex- 
posure time defined in section [3l we study in this section the 
whole telescope transmission i.e. detector sensitivity and filter 
transmission assuming mirror reflectivity of bare silver In sec- 
tion |4T1 we define the filter properties. In sections l4!2l 14.31 and 
14.41 we study the filter resolution to improve the photometric 
redshift accuracy and decrease the number of catastrophic red- 
shifts. 

To develop a survey strategy we need to define the number of 
filters and their shapes, since this will aff'ect the survey speed in 
terms of the required exposure time per filter. In this paper, we 
use square shaped filters and vary the filter shape in changing 
their width. 
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4.1. Properties of filter set 

To optimize the design of the filter set for photometric redshift 
quality, we choose logarithmically spaced filters within a given 
wavelength range (Davis et al. 2006). The wavelength range is 
chosen so as to use the full capacity of the detectors. This log- 
spaced repartition of filters mimics the wavelength shift-dilation 
of galaxy spectra as a function of redshift z, expressed by the 
formula Aghs - Aresti^ + z). The useful spectral features for the 
photometric redshift are shift-dilated as a function of redshift 
and the filter set is designed to follow this evolution allowing a 
direct comparison of galaxy luminosity as a function of redshift. 
Thus, each filter is a redshifted copy of the previous one and its 
width is multiplied by a factor of a (explained below). This filter 
design was first developed fo r the SN probe to improve the K- 
correction jDavis et al. i Booeh . However, this is also relevant for 
photometric redshifts and galaxy evolution studies. We construct 
the first filter using the detector cut-off" in the near-UV of the 
visible CCDs (^"" ~ 3200A) and a width w: (r' , A"' + w). 

• mm ' ^ mm^ mm ' 

The subsequent filters are based on the first filter multi- 
pUed by the factor a. This replication factor is defined as a 
function of the wavelength range available to the instrument 

(i*.f'-,4r'"'0 = = (3200A, 17000A), the 

width of the first filter wo, and the number of filters n 



Filter resolution studies 



^.detector _ n lO 



= ""(KT"" + wo) 



^detector 



\detector 



(26) 



where A^l^^^ is the maximum wavelength of the / = filter 
Following these properties of filter sets, we define the resolution 
of a filter or a filter set as 



A', 



{i=filter) 



^A' 



a'A^ 

mean _ fn 

a' Wo 



■{filter set)- 



(27) 



Using this definition of filter sets, we attempt in Section 14.21 to 
find a filter set resolution that gives the best results in term of 
photometric redshift quality for WL studies using the telescope 
design and noise prescription that we defined in Section 123] 

4.2. Filter resolution studies: methods and hypothesis 



The filter set properties defined in Section HTTI are determined by 
the width wq of the first filter and the total number of filters n. 
This also defines the resolution of the filter set. 

In this section, we study the impact of the resolution on 
the photometric redshift quality and WL analysis. Following the 
studies in Bemtez et al. (2009), we choose to use eight filters that 
they proved to be the optimal number of filters to reach the high- 
est completeness in depth and quality of the photometric red- 
shift distribution. Their studies are based on mock catalogs de- 
rived from the HDF catalog artificially extended to 5000 objects. 
Their phot ometric redshi ft distribution was computed using the 
BPZ code (lBemtej|2000l) . We also study the minimum number 
of filters required using our optimal filter resolution in Section 

m 

To test the impact of the filter resolution, we made a grid 
in resolution by raising wo - the width of the first filter - in 
covering steps of lOOA. Using 16 configurations of filter set 
Wo = [600A, 2000A], we study the evolution of the scatter and 
the number of catastrophic redshifts as a function of the filter 
resolution. 
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Fig. 13. Transmission as a function of wavelength in in the ex- 
fi-eme cases of the filter sets tested: 1R = 6 (top panel) and H-l 
(bottom panel). We include the transmission of optics (long- 
dashed pink), CCD (small-dashed black), NIR detector (dotted 
black), and the four-mirror reflectivities (solid brown). 



Figure [13] shows the two extreme cases of filter resolution 
we tested. The upper panel contains our results for a filter set 
with the highest filter resolution, 1^-6. This high filter reso- 
lution makes the filters very narrow: indeed gaps appear in the 
wavelength coverage, which is something we wish to avoid. We 
note that = 6 is the only filter resolution studied here that has 
wavelength gaps in its transmission. This is a filter configura- 
tion that is not desired unless complemented with broader band 
observations. The bottom panel shows the lowest filter resolu- 
tion H-l with extremely broad filters. The multicolor lines are 
the filter transmission curves. The dashed and dotted lines are, 
respectively, the assumed CCD and NIR detector transmission 
curves. The dot-dashed line is the four-mirror reflectivity curve 
using ba re silver re flectivity as described for the SNAP/JDEM 
mission (iLevill20()7l) . We assume a survey configuration of four 
bare silver mirrors to focus the light on the focal plane using 
a quantum efficiency (QE) similar to the LBNL detector trans- 
mission properties shown Figure [13] These characteristics have 
an impact on the photometric redshift accuracy that depends on 
the photometric errors calculated using equations in section [23] 
For each filter set created, we compute photometric redshifts us- 
ing the Le Phare photometric redshift code brieffy described in 
section 12.11 



4.3. Filter resolution studies: Photometric redshift quality 

Figure [14] shows the photometric redshift scatter cr{\zp - zs\) 
as a function of filter set resolution, binned by magnitude in the 
upper panel and by redshift in the bottom panel. Each square 
point of these curves shows the result for a particular filter set 
configuration with a resolution 2 < !?? < 6. We use the I-band like 
F4 filter for the magnitude binning (see Table [TJ. To minimize 
the photoz scatter a filter resolution of > 3 is preferred when 
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Fig. 14. Photometric redshift scatter cr^ore as a function of filter 
set resolution binned by magnitude up to z 6 (top panel) and by 
redshift (bottom panel) up to I^b - 26 mag. F\ represents the 
I-band filter and Zs the spectroscopic redshift. 



looking at the top panel of Figure[T4] The bottom panel of Figure 
[14] suggests a preferred filter resolution of » 3 - 4. 

The accuracy of the photometric redshifts depends on the 
color gradients of galaxy templates. It also depends on the pho- 
tometric errors that are used as a weight in the template fitting 
procedure (see Equation[T]i. A high filter resolution (% > 5) low- 
ers the S/N in each filter and the weight derived from it do not 
place sufficient constraints to ensure an accurate photoz estima- 
tion. In the case of a low filter resolution (K < 3), the overlap 
between filters lowers the galaxy color gradient degrading the 



quality of the photoz results. Figure [T4l shows that an optimal 
filter resolution is around :R a; 3 - 4. 



4.4. Filter resolution studies: Catastrophic redshift rate 
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Fig. 15. Percentage of catastrophic redshifts as function of the 
filter set resolution binned by magnitude up to z 6 (top panel) 
and by redshift (bottom panel) up to Iab = 26 mag. F4 represents 
the I-band filter and Zs the spectroscopic redshift. 



Figure [H] shows the percentage of catastrophic redshifts as 
a function of the filter set resolution binned by magnitude (top 
panel) and by redshift (bottom panel). Each square point cor- 
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responds to a filter set configuration. Catastrophic redshifts are 
defined as \zphoi - Zspeaml > 0.3. 

To constrain the DE parameters, one of the WL techniques 
consists of dividing the ga l axy distribution into redshift slices 
dBernstein & HuteredllOlOl ISun et al.ll2009 '). We thus need ac- 
curate photometric redshifts to avoid a contamination between 
slices. For the redshift range 1 < z < 3, the Balmer break or 
the D4000 is in the wavelength range fully covered by the fil- 
ter set 3200A < A < 17000A. The color gradient produced will 
thus ensure a robust photoz estimation. In the 2 < z < 3 red- 
shift range, the galaxies are fainter increasing the probability of 
color confusion and resulting in a higher catastrophic redshift 
rate. Consequently, in this redshift range, a higher filter resolu- 
tion increases the color gradient accuracy which improves the 
photoz accuracy as shown by the blue and violet curves in the 
bottom panel of Figure [TSl 

High redshift galaxies (3 < zs < 6) usually have faint appar- 
ent magnitudes. Broader filters are then more suitable to max- 
imize the S/N as shown in terms of the percentage of catas- 
trophic redshifts binned by magnitude (top panel of Figure [TsTl. 
The top panel of Figure[T5]shows a preferred resolution range of 
3 < < 4 with a significant decrease in the fraction of outliers 
for galaxies F4 > 24. It reduces the outlier rate by 3% to 10% 
depending on the magnitude range considered. 

For future dark energy surveys, the WL analysis is based on 
the statistics of faint and numerous galaxies. Figures [14] and [15] 
show that the optimal WL choice uses broad filters and has a 
resolution in the range - 3 - 4. 
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Fig. 16. Number of galaxies and median(|z/7-zi|) as a function of 
the filter set resolution for different photometric redshift quality 
selections. 

Figure [16] shows the number of galaxies and the 
median[z/? - zs] for different photometric redshift quality selec- 
tions as a function of the filter resolution. On the one hand, if a 
strict photoz quality selection is used, the resolution giving the 
highest number of galaxies is in the range of K - 3 - 4. On the 
other hand, the bias is smaller at higher filter resolution. This 
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Fig. 17. Lensing FoM as a function of the filter set resolution 
without Galactic absorption included, and a constant Zodiacal 
light. The dashed blue, dotted red, and solid black curves repre- 
sent respectively a survey area of 10000, 7000, and 3000 deg^. 



shows the importance of an accurate spectroscopic redshift cal- 
ibration to estimate and correct for the bias of the photometric 
redshift distribution. 

To reach a definite conclusion about the filter resolution 
question. Figure [17] shows the FoM (defined in section [3]l as 
a function of the filter resolution for three different fractions of 
the sky observed. A filter resolution of 3.2 provides the tightest 
DE constraints. This FoM calculation does not take into account 
the catastrophic redshift rate. However, this filter resolution cor- 
responds to the lowest rate of catastrophic redshifts (as shown in 
Figure [TSll. hence should be the optimal filter resolution in terms 
of the photoz accuracy of a WL analysis. 

5. CCD's blue sensitivity for catastrophic redshift 

A possible way of reducing the number of catastrophic redshifts 
at low and high redshift is to optimize the efficiency of visible 
detectors in the near-UV. A higher sensitivity in the wavelength 
range [3000-4000A] improves the photoz results at low redshift 
derived from either the Balmer or D4000 break. This results 
in more accurate color gradient and photometric redshifts. In a 
similar way, it also helps to decrease the catastrophic redshifts 
rate at high redshift related to the Lyman break feature. The 
Lyman break is at 912A rest-frame and enters into the filter set 
for galaxies at z ~ 2.5, which will help us to break the color 
degeneracy between low and high redshift galaxies. 

Thus, we explore the impact of the CCD quantum effi- 
ciency (QE) and produce five QE curves differing in the near- 
UV wavelength range. Each detector curve is then used with 
the eight-filter set of a resolution {R) x 3.2 to derive noise 
properties that are applied to generate a realistic mock catalog 
following the prescription described in Section 12.31 using the 
CMC. Photometric redshifts are thereafter calculated using the 
Le Phare photometric redshift code. 
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The bottom panel of Figure [18] shows the five CCD QE 
curves used (solid multicolor lines) as a function of wavelength. 
Each curve defines a mock catalog called "ccdQE.it", where 
X - [0, ..,4]. The blue "ccdQEO" is the most efficient at blue 
wavelengths, having a 40% efficiency at 3620A. The one with 
the worst efficiency "ccdQE4" is shown in purple and has a 40% 
efficiency at 4420A. The numbers below the figure are the wave- 
length at which the CCD QE curves reach 40% efficiency. The 
dot-dashed line is the four mirror reflectivity. We note that the 
four mirror reflectivity produces a cut-off of the bluest detector 
curve. The dashed lines are the four first filters of the eight filter 
set. The first two filters are the most affected by this gain in the 
CCD QE in the near-UV/visible wavelength range. 

Figure[T8]top panel shows the percentage of catastrophic red- 
shifts as a function of the visible detector efficiency at 3600A. 
Each square point is the same eight-filter set convolved with 
a different efficiency of visible detector and each curves rep- 
resent a different magnitude bin. We are mostly interested in 
the magnitude range 25 < < 26.5 since this contains the 
very faint galaxies that will be used in a WL analysis. For 
25 < F4 < 26, the green curve shows that the mock catalog 
ccdQEl contains 13% of catastrophic redshifts, while the cat- 
alog ccdQE4 contains around 21%. A blue-optimised detector 
allows us to strengthen the Balmer break signal in the bluest 
bands, which helps us to decrease the percentage of catastrophic 
redshifts at low-z i.e. < z < 1 and high-z i.e. z > 2.5. 

The ccdQE4 catalog provides simular results to a seven-filter 
configuration covering a wavelength range of [4100 - 17000A]. 
This illustrates how the removal of the first filter would affect 
the percentage of contamination by the catastrophic redshifts in 
a given magnitude range. The fraction of catastrophic redshifts 
is about two times higher in the magnitude range of interest for 
WL, 25 < F4 < 26, when not including the bluest filter. 

We note that there is not much improvement between us- 
ing the catalogs ccdQEO and ccdQEl, which is explained by the 
cut-off in the four mirror reflectivity at 3200A. We note, how- 
ever, that it is critical to improve the blue sensitivity of visible 
detectors (and possibly the reflectivity of the mirrors at blue-UV 
wavelengths) to help remove catastrophic redshifts. Ultimately, 
this can be done by using a detector dedicated to the U-band, 
which would be blue-optimised. 

Figure [19] shows the Zp - Zs distribution of ccdQEO for left 
figures and ccdQE4 for right figures. We computed the photoz 
using a library of SED templates with emission lines for the bot- 
tom panels and without emission lines for the top panels. This 
last configuration is equivalent to a poor photometric redshift 
calibration. The mean and scatter are calculated for the galax- 
ies whose redshift is located inside the two red lines as defined 
in section 12.41 The catastrophic redshift rate (gold writting) is 
the percentage of galaxies outside the red lines. We also calcu- 
lated the scatter cr,,.,,,, for galaxies meeting the 68% confidence 
interval criterion A*^*z < 0.5 as defined in section lZ?] A higher 
efficiency in the 3000 - 4000A range helps to reduce the pho- 
toz scatter and minimize the number of catastrophic redshifts as 
shown in Figure [19] The catastrophic redshift rate is multiplied 
by a factor of two when the U-band photometry is of poor qual- 
ity. In the case of both a poor U-band photometry and a poor 
photometric redshift calibration, the catastrophic redshift rate is 
multiplied by a factor of three as shown in Figure [19] 

To summarize, improving the U-band photometry helps to 
similarly improve photometric redshift estimates and reduce the 
catastrophic redshift fraction by a factor of two in the magnitude 
range of interest for WL. 
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Fig. 18. Bottom panel represents the transmission of optical dec- 
tectors (solid lines), filters (dahsed lines), and mirror reflectiv- 
ities (dot-dashed lines) as a function of wavelength in A. The 
inscriptions below the bottom figure are the wavelength at 40% 
efficiency for the different optical detector curves named ccdQEx 
where x = [0, .., 4]. Top panel represents the percentage of catas- 
trophic redshifts as a function of the efficiency of the optical de- 
tector at 3600A. The colors of points are corresponding to the 
optical detector of the same color that have been used in the 
photometric noise calculation. 



6. WL survey strategy: Number of filters and survey 
area 

Future DE surveys plan to cover large areas to detect a large 
number of galaxies. However, one should carefully consider the 
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Fig. 19. zp-zs distribution in two cases of visible detector effi- 
ciency curves: ccdQEO for panels on the left and ccdQE4 for 
panels on the right. Bottom panels include the effect of emission 
lines in the library of templates used for the photometric redshift 
estimate, while top panels do not. 



survey strategy and the instrument design together to optimize 
the areal coverage. Consequently, one has to choose the pixel 
scale and the number of pixels that determine the FOV of the 
camera, the observation time, the number of exposures, and the 
number of filters. Each of these choices affects the WL analysis 
in terms of the number density of galaxy sources, the photomet- 
ric redshift accuracy, and the shape measurement quality, which 
defines the appropriate number of galaxies. In section [3731 we 
have defined a minimal observation time of 200s for each fil- 
ter. In section l4!2l using this exposure time we found an optimal 
filter resolution of » 3.2 for an eight filter configuration. In 
Figure |20l we show the WL dark energy FoM normalised by 
their respective maximum for a six, seven, and eight filter con- 
figuration as a function of filter resolution. The resolution of 3.2 
provides the highest values of FoM independently of the number 
of filters. We thus use this optimal resolution to compare the effi- 
ciency of six, seven, and eight filter configurations, aiming to de- 
fine a minimum number of filters. With the CMC mock catalogs, 
we simulated three surveys of six, seven, and eight filters using 
a resolution ofR^ 3.2. Figure |2T] shows these filter configura- 
tions. To make a fair comparison between survey configurations, 
we use the same total observation time and divide this by the 
number of filters for a given survey configuration. In all Figures 
of this section, we used a total observation time of 2400^, which 
is distributed into four exposures of respectively 150s, 150s, and 
120s exposures per optical filters for the six, seven, and eight 
filter sets. The 6-filter configuration has four optical filters and 
two NIR, while the 7 -filter and 8 -filter configurations have three 
NIR filters each and respectively four and five optical filters, as 
described in Table |2] The I-band of the three configurations is 
F4 for the 8-filter catalog and F3 for the 7 and 6 filter catalogs. 
However, to make fair selections and look at the same galaxy 
photometric redshift results binned by magnitude, we simulated 
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Fig. 20. Normalised lensing FoM as a function of filter resolution 
for a 6, 7, and 8 filter configurations in respectively dashed blue, 
dotted red, and black solid curves. 



for each survey configuration a magnitude in the F4 filter of the 
8-filter survey configuration and used this band as a reference to 
make different magnitude cuts for all 6, 7, 8-filter survey config- 
urations and allow simple comparisons. 
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Fig. 21. The 6-, 7- and 8-filter configurations using a filter reso- 
lution ofR^ 3.2. We include the transmission of optics (long- 
dashed pink), CCD (small-dashed black), NIR detector (dotted 
black), and the 4-mirror reflectivities (solid brown). 
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The filter set properties are described in the Table |2] where 
!R and Heff are respectively the filter resolution before and af- 
ter convolution with the mirrors reflectivity and detectors QE. 
To have a close effective resolution Hef f for all filters, we con- 
structed the first filter independently of the other filters. We note 
that the detectors QE and mirrors reflectivity have a large impact 
on Heff for the first filter as shown in Table |2] 

Table 2. Filters characteristics. 



Photoz dispersion / mag bins S/N>10 



6.1. Photometric redshift quality vs number of filters 

In Figures I22ll23l and|24l we use the galaxies with S /N > 10 in 
the I-band in the photometric redshift analysis. Figures |22] and 
|23] show, respectively, the photometric redshift dispersion and 
percentage of catastrophic redshifts as defined in section 12.41 
Compared to the 8 -filter configuration, the 7 -filter configuration 
has a photometric redshift scatter that is larger by 0.01, while 
the 6-filter configuration has a dispersion that is 0.02 larger (top 
panel Figure l22ll. However, the percentage of catastrophic red- 
shifts is lower in the 7-filter configuration than in other configu- 
rations as shown in Figure |23] This last figure shows the catas- 
trophic redshift rate as a function of redshift (bottom panel) and 
magnitude (top panel). Both panels show that the 7-filter con- 
figuration is very similar to the 8-filter configuration. However, 
for faint galaxies at magnitude / > 24.5, the 8 filter configu- 
ration has lower catastrophic redshift rates as shown in the top 
panel of figure|23] The 6-filter configuration provides results that 
are relatively worse than other configurations, with about 8% 
of catastrophic redshift contamination for galaxies at magnitude 
/ =s 25.5 ataS/N > 10, while the 7 and 8-filter set have about 
5% contamination in this magnitude range. Even though the 6- 
filter configuration has deeper photometry than other configura- 
tion, this does not compensate for the lack of color gradient in- 
formation (from the amount of filters) needed for a good photoz 

esti mation. 
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does not significantly degrade the estimated DE parameters. 
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Fig. 22. Photometric redshift scatter as a function of the spectro- 
scopic redshift (top panel) and 1-band magnitude (bottom panel) 
for the 6-, 7-, and 8- filter configuration. 



as shown in 'Bridle & Kin^ ( |2007|) . Thus, photoz scatter, bi- 
ases, and catastrophic redshifts may have a strong effect on 
the estimated parameters. Consequently, the recommended 
minimum number of filters is seven assuming our covering 
strategy in wavelength. This configuration gives the highest 
accuracy and minimizes the number of catastrophic redshifts 
binned in redshift and magnitude. We also conclude that the 
6-filter configuration does not seem a good option in terms of 
photometric redshift accuracy. 

Figure |24] shows the cumulative number density of catas- 
trophic redshifts as a function of magnitude (top panel) and red- 
shift (bottom panel). The solid, dotted, and dashed lines show 
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Fig. 23. Percentage of catastrophic redshifts as a function of the 
spectroscopic redshift (top panel) and I-band magnitude (bottom 
panel) for the 6-, 7-, and 8- filter configuration. 



Fig. 24. Density of catastrophic redshifts as a function of the 
spectroscopic redshift (top panel) and 1-band magnitude (bottom 
panel) for the 6-, 7-, and 8- filter configuration. 



respectively the galaxy number density for which i^Zp - Zs\ > 
Q3),{\Zp - Zs\ > 0.3 & A^^'^'z < 0.5), and (A^^'^^z > 0.5). The 7 
and 8-filter catalog have similar results with a number of catas- 
trophic redshifts of 1.1 gal/arcmin^, while the 6-filter configu- 
ration has 2 gal/arcmin^. For the 7 and 8-filter configuration, 
most of the catastrophic redshifts are flagged as a poor photoz 
estimates (A^'^'^'z > 0.5), leaving only 0.3gal/arcmin^ of con- 
tamination when excluding these. For the 6-filter configuration, 
the number density is 0.6gal/arcmin^ after excluding flagged ob- 
jects. Thus, these dotted lines show the number of catastrophic 
redshifts that the code selected as a reliable estimate of galaxy 
redshifts that we define to be A^^'^''z < 0.5 in Section [3] These 



galaxies are the ones that enter into a WL analysis and degrade 
the cosmological parameter estimation. 

This work also demonstrates that for the same amount of 
telescope time the 7 and 8-filter configurations are the best op- 
tions in terms of photoz results. Both of these configurations 
would give an accurate photometric redshift estimate for most 
galaxies down to / a; 26. All these numbers are optimistic since 
we are using the same library of galaxy spectra to generate mock 
catalogs and calculate photometric redshifts. If the photometric 
redshift calibration survey (PZCS) is perfectly representative of 
the galaxy distribution in terms of redshift, depth, color, and size. 
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these results should be very close to what we will obtain in future 
DE misssions. 

6.2. Optimizing tlie WL strategies using tine dark energy FoM 

We use the iCosmo package to compute FoM using the telescope 
characteristics defined in section [3] To calculate the FoM from 
these configurations, we assume a space-based mission with all 
filters onboard the satellite. 

6.2.1. Surveying the sky, limits from Galactic absorption and 
Zodiacal ligint 
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Fig. 25. WL FoM as a function of exposure time by pointing 
with a primary mirror diameter of 1.5m, a pixel scale of 0.15", 
and a FOV=0.5deg^ for a 5-year mission. The magenta curves 
include the effect of the Galactic absorption and Zodiacal light 
variation, while the black curve do not. The dashed curves rep- 
resent the number of galaxiesxle? for both configurations. The 
green region shows the COSMOS incompleteness affecting the 
reliability of CMC galaxy densities, which may impact conclu- 
sions in these areas. 

In Figure |25l we computed the WL FoM (solid lines) and 
total number of objects (dahsed lines) for a five-year mission as 
a function of the exposure time by pointing for an 8-filter config- 
uration. We assume two cameras, with five filters for the visible 
camera and three filters for the IR camera. We assume the same 
total exposure time for the NIR and visible cameras, leading to 
different exposure times for each filter depending on the camera 
the filter belongs to (see Table |2l). The magenta curves illustrate 
the effect of the Galactic absorption and Zodiacal light variations 
as discussed in section [34l eq. l24l while black curves do not. The 
green dashed area shows the CMC incompleteness region, which 
corresponds to the limiting depth of the COSMOS data. Because 
the incompleteness region only occurs at the longest exposure 
time considered, it does not affect the results presented here. This 
green-dashed part of the graph will be affected by the COSMOS 
limiting depth. The axis on the top represents the area in square 



degrees for a five years WL mission, including (magenta lines) 
or excluding (black lines) the variation in both Galactic absorp- 
tion and Zodiacal light. The black curves are what can be found 
if you neglect the Galactic absorption variation across the sky 
(as in lAmara & Refregie^ (12007 ')). The magenta curves repre- 
sent a specific survey strategy consisting of compensating for 
the Galactic absorption and Zodiacal light variation by varying 
the exposure time in a way that the photometry depth is equal 
for the whole area surveyed. This survey strategy has a strong 
impact on the survey speed of the mission. When properly taken 
into account the Galactic absorption and the Zodiacal light vari- 
ation, the often quoted 20,000 deg^ is no longer optimal. The op- 
timal WL survey we found covers ~ 11 OOOdeg^ for a five-year 
mission. 

The survey strategy chosen here is to obtain a similar pho- 
tometric depth across the survey area, thus increasing the expo- 
sure time at low Galactic latitudes. Another possible strategy is 
to use the same exposure time for every pointing, but this would 
result in an inhomogeneus photometric depth across the survey 
area. The photometric redshift quality would then vary across 
the survey area. This would be difficult to compensate for and 
would likely degrade the DE constraints from the WL analy- 
sis as the survey beyond 10000 deg^ will not provide a much 
higher galaxies density than can be achieved from the ground. 
Ultimately, it may be interesting to consider a trade-off between 
survey area and depth in more complexe suvey, especially if we 
include other probes in the FoM calculation, such as BAO or SN. 

We note that the impact of the Galaxy on the survey remains 
simplified, as we did not take into account the variation in the 
star densities. This will reduce the effective surface density of 
extragalactic sky observed at low Galactic latitudes, and support 
the case for an even smaller survey area than the one suggested 
here. 

6.2.2. Photometric redshift quality e 

Figure |26] shows the galaxy number density as a function of 
the quality criterion we use for the photometric redshift selec- 
tion: e - A^^'''''z/(1 -I- zp) (top panel). In the bottom panel of 
Figure |26l we compute the dispersion cr, mean, and median of 
the {zp - zs)l{ \ + zs) distribution as a function of the photomet- 
ric redshift quality criterion. This shows that the 6-filter config- 
uration has a dispersion that is larger by 0.02 to 0.05 than the 
values for the 7 and 8-filter configurations and is also more bi- 
ased as shown by of mean and median of the photometric red- 
shift versus the spectroscopic redshift distribution. We note that, 
again, the 7-filter configuration has similar results to the 8 filter- 
configuration. 

We computed in Figure [27] the WL FoM for a three-year 
mission: 1) as a function of the photometric redshift quality e, 
and 2) as a function of the exposure time by pointing as defined 
in section [J!4l (eq.l24li. and 3) for different numbers of filters. 

The axis at the top represents the area in square degrees for 
a three year WL mission, including the effects of Galactic ab- 
sorption. The magenta, cyan, and blue curves represent the WL 
FoM for respectively the 8-, 7-, and 6-filter configurations as de- 
fined in section 16.11 The solid and dotted lines present the WL 
FOM for respectively a photometric redshift quality selection of 
e = 0.05 and e = 0.1. The green dashed area shows the CMC 
incompleteness region. The incompleteness region varies for the 
different survey configurations. Hence, we present the region for 
the most restrictive case, which corresponds to the 6-filter con- 
figuration in this figure. 
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Fig. 26. Statistics and number density for the 6-, 7-, and 8-filter 
survey configurations using the telescope properties described in 
Table lA.ll with observation times described in Sect. |2] 
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Fig. 27. WL FoM as a function of exposure time by pointing for 
the 6-, 7-, and 8-filter configurations using a primary mirror di- 
ameter of 1.5m and pixel scale of 0.15", and a FOV=0.5deg^ for 
a 3-year mission. The solid and dotted curves represent a photo- 
metric redshift selection of respectively (zph-zpl)/(l-i-zp)<0.05 
and 0.01. 



The solid and dotted curves show that the eight-filter con- 
figuration provides the tighest DE constraints for an exposure 
time by pointing of 3600s and 4000s, respectively. This corre- 
sponds to an exposure time by filters of 180s-200s per filter for 
the visible camera over an area close to 7300deg^ for a three 
year mission. The number density of useful galaxies reaches 45 



gal/arcmin^ at lOcr in I-band for the eight and six-filter con- 
figurations as shown in Figure |28] This last figure represents 
the galaxy number density as a function of exposure time by 
pointing for the three configurations. The six-filter configura- 
tion has results close to the eight-filter configuration, especially 
if we relax the photo-z quality selection. However, this calcu- 
lation does not include catastrophic redshifts, the discrepancies 
between both configurations will be a larger, especially if we in- 
clude the impact of a PZCS and the catastrophic redshift rate. 

The eight-filter configuration yields the best FoM (although 
not by a large factor). The accuracy needed to estimate the DE 
parameters requires an ambitious survey. A 6-filter configuration 
may render the final measurement unreliable because of the sys- 
tematic errors introduced by the larger number of catastrophic 
redshifts and a larger photometric redshift scatter and biases. 
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Fig. 28. Density of galaxies in an I-band filter as a function of 
exposure time by pointing for S/N> 1 (dashed), S/N> 1 & (zph- 
zpl)/(l-i-zp)<0.1 (solid), S/N>10 & (zph-zpl)/(l+zp)<0.1 (solid) 
& galaxy size>1.25xePSF (dotted). 



Figure |28] also shows the impact of the photometric and size 
selections on the three configurations. The most stringent se- 
lection is the quality cut based on galaxy sizes with a loss of 
respectively 8 gal/arcmin^ for the 7- and 8-filter configurations 
and 13 gal/arcmin^ for the 6-filter one. The 6-filter configuration 
has deeper photometry, which raises the number of faint, small 
galaxies thus causing this configuration to be more affected by 
the galaxy size criterion. For the same reason, the 6-filter con- 
figuration is more affected by the photometric redshift selection, 
S /N > 10, and loses around 8 gal/arcmin^, while other configu- 
rations lose around 5 gal/arcmin^. 

In this comparison, we have used by default a 1.5m primary 
mirror diameter If were to consider a 1 .2m mirror diameter the 
loss produced would be a lot more significant because of the S/N 
decrease and the larger galaxy-size selection. 
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Fig. 29. WL FoM as a function of exposure time by pointing for 
a mission time of 1 year (magenta), 3 years (cyan), and 5 years 
(blue). The dashed curves represent the number of galaxiesx le7 
for the three mission durations. 

6.2.3. Effect of the mission duration r„„j„o« 

We explore the influence of the mission duration for a 1 .5m tele- 
scope diameter and an eight-filter configuration with a photo-z 
quality selection of e - 0.1. In Figure l29l we computed as a 
function of the exposure time per pointing the WL FoM (solid 
lines) and the total number of galaxies (dashed lines) for a 1, 3, 
and 5-year mission. The green dashed area shows the CMC in- 
completeness zone. This green-dashed part of the graph will be 
affected by the COSMOS limiting depth. We find that the op- 
timal configuration for a one year survey is a 3000s observing 
time by pointing, which represents 150s per visible filter for a 
total survey area of 3200deg^. For a 3-year mission, the optimal 
observation time per pointing is 3600s, or 180s per visible fil- 
ter, for a total survey area of 7300deg^. For a 5-year mission, 
the optimal observation time per pointing is then 3600s or 180s 
per visible filter for a total survey area of llOOOdeg^. This op- 
timal exposure time corresponds to a number density of useful 
galaxies of 45gal/arcmin^ at lOcr in I-band. The shape of the 
curves indicate that longer missions reach their maximum FoM 
for slightly deeper exposures rather than shallower exposures. 
A one-year mission is likely limited by the RON regime (as the 
Galactic absorption is less important on scale of 3000 square de- 
grees). As the mission duration increases, the Galactic absorp- 
tion becomes more dominant and a deeper photometry provides 
more galaxies than a larger survey area. 

We note that one can trade mission time with the size of the 
field-of-view (number of detectors assuming a fixed pixel scale) 
as discussed in section [3741 For example, a 5-year mission with 
a 0.3 square degrees FOV is in principle similar to a 3-year mis- 
sion with a 0.5 square degrees FOV. 

6.2.4. Impact of tine pixel scale p'" 

Next we compute the impact of the different pixel scales on the 
WL FoM as a function of the exposure time per pointing, assum- 



ing: a 1.5m primary mirror diameter, an 8-filter configuration, a 
5-year mission, and a FOV of 0.5deg^. Figure l30l shows that a 
larger pixel scale results in a larger FoM. We note that the max- 
imum of the FoM is obtained at very different exposure times 
per pointing: smaller pixel scales require longer exposure times 
to reach the photon noise regime. The gain going from small 
to large pixels is significant in terms of FoM, although the gain 
from 0.15" to 0.2" is only 10% in FoM. As there are concerns 
about how well the shape can be measured if the sub-dither is 
not optimal, we assume a pixel scale of 0.15" to be an optimal 
choice. 
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Fig. 30. WL FoM as a function of exposure time by pointing 
for a pixel scale of 0.2"(magenta), 0.15"(cyan), 0.1"(blue), and 
0.05"(green) using an 8-filter set with a resolution of K - 3.2. 
The dashed curves represent the number of galaxiesx le7 for all 
configurations. 



6.2.5. Impact of tine primary mirror diameter Di 

We then consider the impact of the primary mirror diameter on 
the WL FoM results in Figure|3T] We use the optimal pixel scale 
of 0. 15" for a primary mirror diameter of 1 .5m as argued in sec- 
tion [321 and adapt the pixel scale to have the same sampling of 
the ePSF in the 1.2m and 1.8m case. 

An interesting and not intuitive result is that when using a 
primary mirror diameter of 1 .5m, the mission time is better em- 
ployed in going deeper on smaller areas. In 5 -year survey, we 
should then cover about ~ lOOOOdeg^ with 4000sec by pointing, 
which represents 800s by filters in four exposures of 200s (at the 
Galactic pole) for the visible camera. With a primary mirror di- 
ameter of 1.8m, one can cover a wider area (~ 12000deg^) since 
a larger mirror allows us to reach a higher galaxy density in a 
smaller exposure time. With a smaller primary mirror diameter 
of 1.2m, the optimal survey covers a smaller area ~ 8000deg^ 
with longer exposure time per pointing. 

We note that we keep the FoV at 0.5deg^ meaning that the 
number of pixels is not fixed for different telescope diameters. 
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Fig. 31. WL FoM as a function of exposure time by pointing 
for a primary mirror diameter of 1.2m, 1.5m, and 1.8m with a 
pixel scale varying in order to have the same optimal effective 
PSF sampling using a 8-filter set with a resolution ofR- li.l. 
The FOV is fixed at O.Sdeg^, which implies that the number of 
pixels decrease when the pixel scale is larger. The dashed curves 
represent the number of galaxiesxleV for all configurations. 

7. Conclusion 

On the basis of our realistic mock catalogs representative of faint 
galaxies in the Universe (Paper 1), we have optimized the obser- 
vational strategy and instrument parameters of possible future 
DE WL space missions, focussing in particular on the impact of 
the photometric redshift quality on the DE constraints. 

In section 2, we have described the photometric redshift 
technique and photometric uncertainties of a space survey. In 
particular, we have defined the parameters used to quantify the 
quality of the photometric redshifts used in the WL analysis. 

To address the observational strategy of a WL space mis- 
sion, we then defined, in section 3, the key properties of the 
instrument, which consists of a visible camera and an infrared 
camera. We computed the effective PSF, the overall throughput, 
and the sensitivities, which depend on the detector characteris- 
tics, the pixel size, the exposure time, the dithering strategy, and 
the value of the sky background. In particular, we investigate the 
impact of the pixel scale on the sampling of the effective PSF 
and place upper limits on the pixel scale values. We also sug- 
gest an optimal pixel scale based on simple arguments, which is 
further investigated in section 6 of the paper. 

In section 4, we conducted a detailed analysis of the opti- 
mal filter set to minimize the photometric redshift scatter and bi- 
ases along with the number of catastrophic redshifts. We found 
that square filters of resolution ~ 3.2 maximize the photomet- 
ric redshift quality for the faint galaxy population used for WL 
measurements in the cases of 6, 7 or 8-filter sets. For a resolu- 
tion of'R~ 3.2, the filter width maximizes the S/N along with 
the galaxy colors resulting in an optimal photometric redshift 
quality. We showed that this H ~ 3.2 resolution maximizes the 
FoM of the WL tomography irrespective of the survey area of 
the mission. 



In section 5, we addressed the issue of the blue sensitivity 
of the visible camera by looking at the quality of the photomet- 
ric redshift and in particular the catastrophic redshift rate. We 
found that improving the blue sensitivity can reduce the catas- 
trophic redshift rate by a factor of two in the magnitude range 
we are interested in for WL analysis. The detector blue sensitiv- 
ity is then a critical parameter since it will help reduce biases in 
the WL cosmological constraints coming from the photometric 
redshift quality. We propose a minimum throughput of 30% at 
3600 A for the visible camera. 

Finally, in section 6, we investigated the impact of some key 
instrument parameters and the observational strategy to maxi- 
mize the FoM of a WL survey. In particular, we showed that 
the instrument parameters and survey strategy are both impor- 
tant in optimizing a WL space mission and that they cannot be 
addressed separately. We thus defined the survey strategy com- 
puting the survey area as a function of the exposure time per 
pointing, the field-of-view of the instrument, the mission dura- 
tion, and the sky model, which includes both the Galactic ab- 
sorption and Zodiacal light variation. 

Quantitatively, assuming a 1.5m primary mirror telescope 
and investigating the ePSF sampling, we argued that a pixel scale 
of 0.15" for the visible camera is the optimal choice to conduct 
an efficient WL survey. This choice is the best trade-off value 
that is small enough to ensure a sufficiently well-sampled ePSF 
after dithering and maximize the FoM unlike smaller pixel size. 
However, this result is based on simple arguments and assumes 
that the PSF of the instrument is sufficiently stable and can be 
well calibrated - this should be a strong requirement of the mis- 
sion and needs an in-depth optimisation. 

Furthermore, assuming a 0.5 square degree field-of-view for 
both visible and infrared cameras, and an eight-filter configura- 
tion with 2l R - 3.2 filter-set, we demonstrated that a ~11 000 
sq.deg survey reaching a homogeneous depth of I/ib=25.6 (at 
lOcr) (which can be achieved with four exposures of ~200s per 
filter, at the Galactic poles) is the optimal and safest combina- 
tion to maximize the weak lensing FoM. The survey strategy 
consists of varying the exposure time to keep an equal galaxy 
number density across the survey area. At this depth, the effec- 
tive number density of galaxies that can then be used for WL is 
45gal/arcmin^, a factor of two better than a WL ground-based 
survey. 

In particular, we show that: 

- An 8-filter configuration is better than a either 6 or 7 -filter 
configuration, which both provides poorer photometric red- 
shift quality and lead to a larger number of catastrophic red- 
shifts. 

- The proper calculation of the survey strategy, including the 
Galactic absorption and the Zodiacal light variations as a 
function of the position of the sky, strongly limits the opti- 
mal survey size of the WL probe to a maximum of ~ 1 1 000 
sq.deg. The exact optimal survey size will depend on the to- 
tal mission time, the pixel scale of the visible detector, the 
primary mirror diameter of the telescope, and the exposure 
time per pointing, which should all be adapted accordingly. 

These conclusions are drawn from a survey strategy that en- 
sures an homogeneus depth across the survey area. Different 
survey strategies may be interesting to consider. A trade-off be- 
tween survey depth and sky area covered might lead to improved 
dark energy constraints. We stress that our analysis is still rela- 
tively simplified, and that further work is needed to investigate 
some important issues. In particular: 1) our figure of merit re- 
main simple, as they do not take into account errors produced 
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by catastrophic redshift; 2) the pixel size optimization is done 
here using simple arguments, and would benefit from a proper 
analysis conducted in the global optimization scheme described 
in this paper; 3) the high number density of stars at low Galactic 
latitude will impact the useful survey area and thus the number 
of useful galaxies to be used in the WL survey, which may shift 
the optimal area to even smaller size. 

To conclude, our analysis addresses the complex optimiza- 
tion of future WL DE space survey by including both the instru- 
ment parameters and the observational strategy. Since some of 
our results were unexpected, we believe that a full optimization 
that includes both the instrument parameters and the observation 
strategy is required to maximize the cosmological constraints of 
the future DE space mission. 
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Appendix A: Typical noise properties for a space 
based survey 

A.1. The galaxy flux e,ig„ai 

To compute the galaxy flux esig„ai, we derive the total flux re- 
ceived by the telescope and the distribution of the light over the 
detectors. The former is given for each band by the magnitude in 
the mock catalog, while the latter depends on the galaxy profile 
and the observed effective radius r^ff, as well as the resolving 
power of the telescope. We assume for all galaxies an intrinsic 
exponential profile following: 

I(r) = exp{^ - L6785.(^)j. (A.l) 

In analyzing the Sersic indices (Sersic 1968) of the galaxy pro- 
files of the zCOSMOS survey dTasca e t al. 2009), we found that 
most galaxy profiles tend to follow an exponential profile (Sersic 
index close to n=l) as shown in Figure IATT] This is particularly 
true for the faint galaxies that represent the bulk of the galaxies 
used for WL measurement. 

The resolving power of the telescope is the size the Airy disk 
of the telescope and its FWHM is defined by 

FWHM[Airy Disk] = FWHM,ei = 1.02.—, (A.2) 

where A is the mean wavelength in the filter considered and Di 
the diameter of the primary mirror. We assume that the FWHM 
of the galaxy can be roughly derived from the convolution of five 
Gaussians representing: 1) the PSF of the telescope FWHM,ei, 
2) the pixel scale p = pir/ccd^ 2) the pixel diffusion cr^, 4) the 
jitter of the telescope cry (see Table lA.ll . and 5) the galaxy theo- 
retical FWHM Ixreff. Thus the observed FWHM of a galaxy is 
expressed as 

FWHMfj = 
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We then need to derive the number of pixels corresponding to 
fracphoi, which corresponds to the number of pixels within the 
circular area of radius 1 .4 FWHM"'" 

gal 



N, 



pvisjir 

0.1 

(A, 



pix 



1.4 FWHM" 



(A.7) 



We can thus derive the number of electrons from the Zodiacal 
light in a similar way as we did for the number of electrons pro- 
duced by the galaxy flux Csignai 



-sky 



Jsky 



Tobs-^ -^-N expo pix 



(A.8) 



where {Q.,Nexpo,Tobs,S) are defined in Table lA. H and 77 is the to- 
tal transmission of the telescope. We note that the diff'erence be- 
tween the number of electrons produced by both a galaxy esig„ai 
and the Zodiacal light e^ky, lies in the terms fracpkot for galax- 
ies and the term Npix-O. for the Zodiacal light. We assume a total 
galaxy flux after being affected by the instrument efficiency for a 
given pixel scale, PSF convolution, jitter that decreases the flux 
by a fraction of it fracphot, while the total Zodiacal light is com- 
puted in the circular area defined by 1 .4 FWHM°^'J correspond- 
ing to NpixO. arcsec^. 



where the term containing the pixel scale contains a 0.5 factor 
representing the combination of the rms size of a box and the 
2.36 factor needed to convert the rms in FWHM. We defined the 
fraction of photons fracphot as the luminosity enclosed inside 
we calculated using the galaxy intensity pro- 
file defined by Eq. lA.ll The fraction of photons received is then 
converted into the number of electrons from the source. Thus, 
the total number of electrons is 



A.3. The detector noises 

The detector noises correspond to the dark current from the ther- 
2 mal radiation of detectors and the read-out noise due to the elec- 
tron motion when reading the detector. The dark current is thus 
'expressed in e^/pix/s, while the read-out noise is in / pix 
3-)Since this is not time dependent. The visible and IR detectors 
are a compound of different materials and thus have different 
characteristics. The visible detectors are CCDs or charge cou- 
pled device using semi-conductor materials such as silicon and 
work mainly in the visible wavelength range. The HgCdTe are 
IR detectors with mercury cadmium telluride. The wavelength 
response can be varied by adjusting the alloy composition of this 
ternary compound. Both detectors need to be cooled to decrease 
both types of noise. The values of the dark current and read-out 
noise are listed in Table I A. fl 



^signal — y signal- 



'oh. 



.S.T].N,_ 



.frac 



phot 5 



(A.4) 



where Nexpo and S are listed in the Table lA.ll and y signal is the 
total number of galaxy photons arriving at the telescope. Tots 
the observation time, and rj the total transmission including the 
telescope, the mirror reflectivity, the detector QE, and the filters 
and optics transmissions, whose values are given in Table [1] 



A.2. The background flux e^ky 

To calculate the sky emission eskx, we used the Zodiacal light 
parametrisation described in Aldering (2002), which gives a use- 
ful approximation of the sky background fluxes using a broken 
log-linear relation in erg cm^^ s^^ A^^ arcsec^^ 



logw(f(A)) = 17.755 for 0.4 < ^ < 0.6^m, (A.5) 
logm(fW) = 17.755 -0.73(/l-0.61;um) for 0.6 < i < 2.2. 

(A.6) 
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Table A.l. Telescope characteristics studied. 



Quantity 


abbrev(if necessary) 


Values 


Primary mirror 




1.5m 


Secondary mirror 


Di 


0.6m 


Observation time 




200s 


Field of view 


FOV 


0.5deg2 


Nbre of exposure 


N 

^ ' expo 


4 




C 

o 


2 j I 2 j j 


Pixel solid angle 


Q. 




Telescope jitter 


o-j 


0.02" 


Pixel diffusion scale 




0.4 pixel 


Pixel diffusion size 




2.6// 


CCD pixel scale 


p 


U. 1 J 


CCD pixel size 




10.5yum 


CCD read-out noise 


„vis 
''RON 


6 / pix 


CCD dark current noise 


P" 


0.03 e^/pix/s 


HgCdTe pixel scale 


0.26" 


HgCdTe pixel size 




18/zm 


HgCdTe read-out noise 


^RON 


5 e'/pix 


HgCdTe dark current noise 


''dark 


0.05 e'/pix/s 



